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Abstract 

In this paper, we study the output feedback stabiHzation for a scalar conservation 
law with a nonlocal velocity, that models a highly re-entrant manufacturing system as 
encountered in semi-conductor production. By spectral analysis, we obtain a complete 
result on the exponential stabilization for the linearized control system. Moreover, by 
using a Lyapunov function approach, we also prove the exponential stabilization results 
for the nonlinear control system in certain cases. 
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1 Introduction 

In this paper, we study the scalar conservation law 

pt{t,x) + {p{t,x)X{W{t)))^=0, tG (0,+oo),xe (0,1), (1.1) 

where ^ 

W{t) = / p{t,x)dx. (1.2) 
Jo 

We assume that the velocity function A is in C'^(M; (0,-|-oo)). Let us recall that the special 
case 

X{s) = -—, sG[0,+cx)), 
1 + s 

was, for example, used in [6} [26]. 

In the manufacture system, the initial data is given as 

p{0,x)= po{x), XG(0,1). (1.3) 
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For this control system, the control is the influx 



uit) := p{t,0)X{W{t)) 



(1.4) 



and the measurement is the outflux 



y{t) ■.= p{t,l)X{W{t)). 



(1.5) 



We consider the following output feedback law 



u{t) - p\{p) = k{y{t) - pX{p)), t G (0, +oo), 



(1.6) 



in which A; € M is a tuning parameter and p G M is the equilibrium that we want to stabilize 
when time t goes to +oo. 

The conservation law that we study here is used to model the semiconductor manufactur- 
ing systems, see e.g. [SI [MIES]- These systems are characterized by their highly re-entrant 
feature with very high volume (number of parts manufactured per unit time) and very large 
number of consecutive production steps as well. The main character of this partial differential 
equation model is described in terms of the velocity function A which is a function of the 
total mass W{t) (the integral of the density p). 

The control problems for conservation laws and general hyperbolic equations/systems 
have been widely studied. The controllability of nonlinear hyperbolic equations (or systems) 
are studied in [lOl O [22l [281 [291 [36] for solutions without shocks, and in [B [21 [3l H [71 
[T ^ \25 \ [52] for solutions with shocks. As for asymptotic stability/stabilization of hyperbolic 
equations (systems), two main strategies have been used. The first one relies on a careful 
analysis of the evolution of the solution along the characteristic curves; see in particular 
O [3 [201 [271 [311 [33] . The second one relies mainly on a Lyapunov function approach; see, in 
particular, [IIl[l2l[l3l[I71[l8l[2ll[35ll37]. 

Concerning the manufacturing model of (jl.l|) itself, an optimal control problem, motivated 
by ini [26], related to the Demand Tracking Problem was studied in [15] (see also [33] for a 
generalized system where A = X{x,W{t))). The objective of that optimal control problem 
is to minimize, by choosing influx u{t) = h{t) instead of ()1.6p . the L^-norm (p > 1) of the 
difference between the actual out-flux y[t) := p(t, 1)X{W {t)) and a given demand forecast 
ydit) over a flxed time period. This is an open-loop control system. Another related work [S], 
which is also motivated in part by [6l|26], addressed well-posedness for systems of hyperbolic 
conservation laws with a nonlocal velocity in . The authors studied the Cauchy problem in 
the whole space without considering any boundary conditions and they gave a necessary 
condition for the possible optimal controls. In a recent paper |16j . controllability of solution 
and out-flux for (|1.1|) have been obtained by the same authors of this one. 

In this paper, we study the exponential stability of the manufacturing system (jl.ip under 
the feedback law (jl.6p . The problem of exponential stabilization can be described as follows: 
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For any given equilibrium p € M and any initial data pQ, does there exist A; € M such that 
p is exponentially stable for the closed- loop control system (jl.ip . (jl.3p together with (|1.6p . 
namely, the weak solution p to the Cauchy problem (jl.ip . ()1.3p and ()1.6p converges to p 
exponentially when time t goes to +00? 

If p = 0, the situation is simple. Natural feedback laws can drive the state to zero 
exponentially fast. In particular, the zero control produces a solution which vanishes after a 
finite time. Nevertheless, if p 7^ 0, the situation is much more complicated. More precisely, 
the stabilization results depend on the equilibrium p 7^ and the velocity function A through 
the following quantity 

• A(p) • ^'-'^ 
First we establish a complete result for the linearized system by spectral approach. A 
sufficient and necessary condition of exponential stability of the linearized system is given 
in Theorem 13.11 However, an example (see [231 page 285]) shows that an arbitrary small 
perturbation of the characteristic speeds of a linear hyperbolic system may break the stability 
property, hence it seems difficult to deduce the exponential stability of the original nonlinear 
problem from the exponential stability of the linearized system. In order to overcome this 
difficulty, we also use Lyapunov function approach to prove exponential stability for the 
linearized system and, then, use the same Lyapunov function to prove the (local) exponential 
stability for the nonlinear system. The Lyapunov functions that we construct in this paper 
are inspired by [21 [TBI [57] . However they have to be modified according to the nonlocal 
feature of the nonlinear system. 

The structure of this paper is as follows: In Section [2l we prove the well-posedness of the 
nonlocal closed-loop system. Then, in Section [3l we study the exponential stability of the 
linearized system. The main results on the stabilization of the nonlinear problem. Theorem 
14.11 and Theorem 14. 2[ and their proofs, are given in Section [H 



2 Preliminaries 

In order to stabilize the system by feedback controls, we first need to recall the usual 
definition of a weak solution to the Cauchy problem (see, e.g., |11[ Section 2.1]) and then 
prove that the closed-loop system is well-posed. 

Definition 2.1. Let p e M, p G [1, +00), /c € M and po € L^{0, 1) be given. A weak solution 
of the Cauchy problem 

pt{t, x) + (p(t, x)X{W{t)))^ = 0, t£ (0, +00), X e (0, 1), 
< p(0,x) = po(x), x€(0,l), (2.1) 
u{t) - pX{p) = k{y{t) - pA(p)), t e (0, +00) 
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+ 



is a function p e C'^([0, +00); i^^(0, 1)) such that, for every T > 0, every r e [0, T] and every 
ip G Ci([0,r] X [0, 1]) such that 

¥?(t,x) = CVx G [0, 1] and (^(i, 1) = 0, Vi G [0, r], 

one has 

1 .1 

x) + A(VF(t))992:(t, x))dx(it — / po(a^)v^(0, x)(ix 
^0 JO 

1) - [ky{t) + (1 - fc)pA(p)](/p(t, 0)) dt = 0. 

Then, the foUowing Lemma holds. 

Lemma 2.1. Let p G M, p G [1, +00) and /c G M 6e given. For any given pQ G -^^(0, 1), i/ie 
Cauchy problem ()2.ip /las a unique weak solution p G C'^([0, +00); LP(0, 1)). Moreover, for 
every T > 0, i/ie foUowing maps 

po GL*'(0,l)^pGC°([0,r];Lf(0,l)), (2.2) 
po G LP(0, 1) ^ {u, y) G L^'IO, T) x LP(0, T), (2.3) 

are continuous. 

Proof. The proof of Lemma l2.1l is similar to that of Theorem 2.3 in [15], thus we only show 
the main ideas and omit the details of the proof. 
Introduce the characteristic curve: 

^ = A(H^(s)), s>0, 

where W{s) = p{s, x)dx. Since p is constant along the characteristics, one can define a 
solution candidate in terms of ^ as following 

'po(i-«t)), itO<{(t)<x<l, 

if < X < i{t) - n+ l<l, orO< i{t) - n < x < 1 for n G N. 

Then by contraction mapping principle and fixed point argument as in |15j . one can prove 
that (|2.4p is indeed the unique weak solution to the original system (j2.1|) . Thanks to (|2.4p . 
one can get the continuity of the maps (|2.2p and (|2.3p : See the proof of [Ml Theorem 4.1]. 

□ 

Our next lemma is straightforward and we omit its proof. 

Lemma 2.2. Let p G M and G M 6e given. If pQ ^ ^'"'^([0) 1]) satisfies the compatibility 
conditions 

A( / po{x)dx){po{Q) - kp^{l)) - (1 - k)-p\{-p) = 0, 

°i 1 
A( / p^{x)dx){p'^{Q) - kp'^{l)) - A'( / p^{x)dx){po{^) - po(l))(/3o(0) - kp^{l)) = 0, 

then the Cauchy problem (j2.ip admits a unique classical solution p G C'^([0, +00) x [0, 1]). 
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3 Stabilization to p for the linearized system 

Before studying the nonlinear control system (j2.ip . we first linearize it near p E M and 
then study the linearized closed loop system: 

Pt{t,x)+\{-p)p^{t,x) =0, t € (0,+oo),x G (0,1), 
< p(0,x) =po(x), xG(0,l), (3.1) 
u{t) = ky{t) + {k- l)dW{t), t G (0, +oo), 

where d is given by (|1.7p and 

W{t):= [ p{t,x)dx, u{t) ■.= p{t,0)X{p), y{t) ■.= p{t,l)X{p). 
Jo 

The stability result for (j3.ip can be stated as follows: 

Theorem 3.1. Lefp G M 6e a constant. Then, G L^{0, 1) is exponentially stable in L^(0, 1) 
for the closed loop system (|3.ip if and only if d > —1 and \k\ < 1. That is to say: if and only 
if d > —1 and \k\ < 1, there exist constants C = C{'p,k) > and a = a('p,k) > such that 
the following holds: For any po G L'^{0, 1), the weak solution p G C^CfO, +oo); L^{0, 1)) to the 
Cauchy problem ()3.ip satisfies 

\\p{t, OIIl^co.i) < Ce"°*IIPo||L2(o,i), yt G [0, +oo). (3.2) 

Next, we prove Theorem 13.11 by spectral analysis in Section [3.11 While in Section [3.21 we 
give another proof, relying on a Lyapunov function approach, that d > —1 and |fc| < 1 imply 
that G L'^{0, 1) is exponentially stable for the closed loop system (j3.ip . The Lyapunov 
functions, which are constructed in Section 13.21 will be used for the stabilization of the 
nonlinear control system as well. 

3.1 Proof of Theorem 13.11 by spectral analysis 

Without loss of generality, we assume 

m = 1- 

Otherwise, a scaling transformation t can easily make it. Thus, by ()1.7p . 

d = pX'{p). 

Then we omit the ~ symbols in Section [3T] and Section [3T2l and rewrite (|3.ip into the following 

Ptit, x) + p^{t, x)=0, t G (0, +c5o), X G (0, 1), 
< piO,x)=poix), XG(0,1), (3.3) 
_p(t, 0) = kp{t, l) + {k- l)dW{t), t G (0, +oo). 
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where W{t) = p{t, x)dx. The exponential decay estimate ()3.2p turns to 

l|p(tr)|lL2(o,i) < C^e-"*||po||L2(o,i), Vt E [0,+oo). (3.4) 

Applying the results in [30], we have the following propositions. 

Proposition 3.1. Let /c € M and pQ € L^(0, 1) 6e given. The Cauchy problem (|3.3|) /las a 
unique solution p € C°([0,+oo);L2(0, 1)). 

Proposition 3.2. Let S{t) {t > 0) be the Cq semigroup on L^(0, 1) that corresponds to the 
solution map of (j3.3p and A be the infinitesimal generator of the semigroup S{t) {t > 0). 
Denote o'p{A) and cr{A) as the point spectrum and the spectrum of A, respectively. Then, 

ap{A) = a{A). 

Proposition 3.3. Let uj{A) := mf{uj G M | 3M = M{uj) : ||S'(t)|| < A/e'^*, Vt > 0} and 
s{A) := sup{9^(;u) I p S a{A)}, where ^{p) denotes the real part of p. One has 

uj{A) = s{A). 

Obviously, by Proposition 13.21 and Proposition 13.31 Theorem 13.11 is equivalent to the 
following Lemma: 

Lemma 3.1. One has s{A) < if and only of d > —1 and \k\ < 1. 

Proof of Lemma l3.lt We only need to study the eigenvalues of the system (|3.3|) . Let G C 
be an eigenvalue of the system ()3.3p and (p ^ he a corresponding eigenfunction. The pair 
{p, (j)) satisfies 

L<^(x) + <^'(x) = 0, xG{0,l), 
|(/)(0) = k^{l) + {k- l)d cp{x)dx. 
For ;U € C, the existence of 7^ such that (|3.ip holds is equivalent to 

1 - ke-^ + (1 - k)d [ e^^dx = 0, (3.5) 
Jo 

the corresponding eigenfunction being (j){x) = e~^^ (up to a multiplicative factor). 
Then we analyze the solution for the characteristic equation (|3.5p in various cases. 
Case 1. d = —1 and A; G M. 

Obviously, ()3.5p admits a zero eigenvalue /i = if and only if 

(1 + d)(l - fc) = 0. (3.6) 

Hence, if d = —1, ()3.5p admits a solution /i = which shows immediately that G -^^^(0, 1) 
is not asymptotically stable for (|3.ip whatever /c G M is. 
Case 2. d ^ —1 and k = 1. 
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In view of ()3.6p . ()3.5p admits a solution fi = and thus € L^(0, 1) is not asymptotically 
stable for 

Case 3. d ^ —1 and k ^ 1. 

We need to analyze the nontrivial solution of the following equation: 
1 - ke-'' + d{l - k) ^ ~ ^ ^ =0, 

It is equivalent to study the zero points of the following continuous function: 

( 1 - 

1 - ke~'' + d(l - k) , if 7^0, 

fdAf^) ■■= \ ^ (3.7) 

{l + d){l-k), if^ = 0. 

Case 3.1. d ^ -1 and \k\ > 1. 

We will apply degree theory for homotopic functions (see Appendix B]) to show that 
fd,k{fJ-) has infinite zero points in the right half plane {fi S C|K(/i) > 0}, and, therefore, 
G i^(0, L) is not stable for (j3.ip . In fact, fd^k behaves close to 1 — ke~'^ as \fi\ — > +oo. 

Let 

H{e, d, k, fi) := fedM = 1 - fce-^ + dd{l - k) ^ ~ \ ^^0. (3.8) 

Then, in particular, H{0, d, k, fi) = /o,fc(/i) = 1 — ke~^ and H{l,d, k, fi) = fd,kifJ')- Obviously, 
/o,fc vanishes at /i^ „ with 

{In/c + i2mT, if k > 0, 
Vn e Z. 
ln|/c| +i(2n + l)7r, if k < 0, 

Note that, since \k\ > 1, ^fik,n > 0. For any fixed n G Z and e > 0, let 

^l,n £ Cjll - A;e~''| < e and |^ - fik,n\ < 1} C C. 

It is easy to see that r2| ^ is a bounded open set of C and 

deg(/o,fc(/i),f^|,„,0) = l, VneZ. (3.9) 

One also easily checks that, if e > is small enough, 

dni^C{n£C\\l-ke-^'\ = e}, Vn e Z, (3.10) 
ni^^ C{fie C\n{n) > 0}, Vn G Z. (3.11) 

We now fix e > small enough so that (|3.10p and ()3.1ip hold. Notice that, for every G [0, 1] 



and any fi G dftf 

i + l±£ 

\H{e,d,k,n)\ = \1 - ke-f" + ed{l - k)- — —\ > e - 1^(1 - A:)||d|^ 
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Hence, for any fixed e > 0, there exists N G such that: for all n € \n\ > N, 

\H{e,d,k,^i)\ > I > 0, yee [o,i],v^ g dn%^^. 

Applying degree theory |1H Appendix B], we get for all n S Z, |n| > iV and all 9 € [0, 1] that 

deg{fdMM,n, 0) = deg{H{l, d, k, 0) = deg{H{e, d, k, /i), n%^^, 0) 

= deg(ii'(0,d,A:,/x),f)|,„,0) = deg(/o,fc(/i),J^|,„,0) = 1. 

Therefore, fd^k has one zero point in r2| ^ C G C\^{p) > 0} for every n G Z, |n| > A'". 
Case 3.2. d < -I and -I <k <l. 

Notice that fd,k{0) = {I + d){l — k) < and that fd,k{p) — > 1 as ^ ^ +oo with ;U € M. 
These facts together with the continuity imply that fd^k has at least one zero point in (0, +oo). 
Case 3.3. d> -I and k = -1. 

We prove that in this case fd-i has infinite zero points on the imaginary axis. Let 
^ = i5 (6 G M \ {0}) be such that 

1 _ p-ib 

fd,-iiib) = l + e-''' + 2d 77— = 0, (3.12) 

lb 

i.e., 

ib{l + cos b — i sin b) + 2d{l — cos b + i sin b) = 0. 
Hence fd-i{ib) = if and only if 

6(l + cos6) + 2(isin6 = 0, (3.13) 
6sin6 + 2d(l -cos6) = 0. (3.14) 

Note that (j3.13p and ()3.14p together with 6/0 imply that cos 6 ^ 1. Then, using also the 
identity 

6(1 + cos 6)(1 — cos 6) = 6sin^ 6, 
one gets that 6 G M \ {0} is a solution of ()3.12p if and only if 6 G M \ {0} satisfies 

6 sin 6 ,66, , , 

5 ^ •= V( — ir-r\ -d = --cot--d = o. 3.15 

2(cos 6 — Ij 2 2 

Obviously, for any fixed d > —1, g{-) G C^{2n'K, 2{n + l)7r) for all n G Z and 

g{b) — )• — oo, as 6 — ?> 2n7r~, Vn G Z"*", 
5(6)^+00, as 6 ^ 2(n + l)7r+,Vn G Z+. 

Therefore, g{-) vanishes at least once in each interval {2mr,2{n + l)7r) for n G Z+, which 
implies that fd-i has infinite zero points on the imaginary axis. 
Case 3.4. d> -1 and \k\ < 1. 

We apply degree theory |1H Appendix B] again for homotopic functions to show firstly 
that s{A) < 0, namely, fd,k has no zero points in the right half plane {/i G C|K(/.f) > 0}. 
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Let H{e, d, k, fi) be defined by ([SSD and for R>0, 

:= {/i G C\^{fi) > and < R}. 

For any i? > 0, H{0, d, k, fi) = /o,fc(^) has no zero points in Qr since |/o,fc(/u)| > 1 — \ke~^\ > 
1 - |A;| > for all ^ G C, > 0. 

Then we claim that: for R > sufficiently large 

H{e,d,k,fi) = fedM / 0, V0 G [0,1], G {/i G C|K(/i) > and |^| > R}, (3.16) 
Property ()3.16p readily follows from 



\fedM\ > 1 - |A;| - 2 ^\ V/i G C \ {0} such that > 0. 

Next we claim that: for R > sufficiently large, 

H{9, d, k, fi) = fedM / 0, V/i G dQR. (3.17) 

To prove ()3.17|) . one first points out that, by ()3.16p . it is sufficient to prove that fed,k does 
not vanish on the imaginary axis. We use a contradiction argument to prove this fact. 
Let /i = z6 (6 G M \ {0}) be a zero point of fed,k, i-e., 

1 — 

Mfc(i5) = 1 - /ce-*^ + 9dil - k) -— = 0. (3.18) 

10 

Property ()3.18p is equivalent to 

5(l-/ccos6) + 6ld(l-/c)sin6 = 0, (3.19) 
- khsiub + ed{l - k){l - cosb) = 0. (3.20) 

Multiplying ()3.19p with 1 — cos h and ()3.20p with sin 6, one gets that 

6(1 - A; cos 6) (1 - cos 6) = -kbsin^h, (3.21) 

Equality ()3.2ip is equivalent to 

6(1 + A;) (1 -cos 6) = 0, 

which implies cos 6 = 1, and thus sin 6 = 0, since 6 7^ and \k\ < 1. Substituting cos 6 = 1 
and sin 6 = into ()3.19p . we get that k = 1. It is a contradiction with the fact that < 1. 
This concludes the proof of ()3.17p . 

Property (|3.17p and the degree theory for homotopic functions for H{9, d, k, fi) give that 

deg(/d,fe(/i), ^R, 0) = deg(i?(l, d, k, n),^^, 0) 

= deg(i/(0, d, k, fi),nR, 0) = degifoM, flR, 0) = 0. 
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Therefore, fd^k dose not vanish in and further in the right half plane {fi € C\^{p) > 0}, 
namely, s{A) < 0. 

Finally, we show that s{A) < 0. For any fixed d > —1 and k G ( — 1, 1), there exists r > 
such that 1 - /ce"^ > 1 - |A;||e-''| = 1 - |/c|e-^(^) > for all e {fi e C\^{fi) > -r}. Since 
fd,kii^) ~ 1 — /ce~'' as tends to +oo, there exists then R > such that fd,kif^) 7^ for all 
p G {fi £ C\^{fi) > — r, 1^1 > R}. If s{A) < —r, we are done. Otherwise, — r < s{A) < 0, 
then s{A) must be achieved by some fi £ {p £ C\fd,k{fJ') = 0, \fi\ < R,—r < ^(/i) < 0} since 
/i I—)- fd,k{p) is continuous. Note that fd,k{p) = has no solution on the imaginary axis, 
we conclude that s{A) < which concludes the proof of Lemma l3.ll and thus of Theorem 

[xn □ 

3.2 Proof of Theorem 13.11 for the case d > —1 and |A;| < 1 by a Lyapunov 
function approach 

Here, using a Lyapunov function approach, we give another proof of the exponential 
stability if d > —1 and \k\ < 1. Note that the solution map of Cauchy problem (|3.3p defines 
a Co semigroup S(t) (t > 0) on L^(0, 1), without loss of generality, it suffices to construct 
the Lyapunov function for every classical (i.e. C^) solutions. An important fact is that, as 
we will see in Section [U the same Lyapunov function also works for the nonlinear closed loop 
system. 

We divide our proof into two cases: |c?| < 1 and d > 1. 
Case 1: \d\<l. 

We construct a Lyapunov function as follows: 



L{t):=[ e-^''p^{t,x)dx + aW^{t), Vi G [0, +oo) 
Jo 



(3.22) 



where the constants /3 > and a G M are chosen later. (The introduction of e is motivated 
by [Hi [131 [351 [37|.) By the definition of W(t) and the Cauchy-Schwarz inequality, we know 

W\t)< Ce'^'^dx C e-'^-'p^{t,x)dx = ^-^ [\-^''p\t,x)dx. (3.23) 
JO Jo P Jo 

If 

a>-^, (3.24) 

L{t) is positive definite for all t > and there exists two constants Ci = Ci{a, f3) > (i = 
1,2,3,4) such that 

Ci||p(*,-)lli2(oi) <C2 j\-^^p\t,x)dx<L{t) 

^ (3.25) 
- ^'i ^ C'4||/0(t,-)|li2(o,i), Vt G [0,+oo). 

Let us first compute the time derivative of L along the solution of (|3.3|) . Note that 

W{t) = - [ p^{t,x)dx = p{t,0)-p{t,l). (3.26) 
Jo 
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It follows from iK22^ . iK2^ and ([33]) that 



L{t) = - / e-'^^(p2(t, x))^dx + 2aW{t)W{t) 
Jo 

= -/3 [ e-^'^p'^it, x)dx - le-^'^p^it, x)] """^ + 2aW{t){p{t, 0) - p{t, 1)) 

= -/3 e-^V^(i, + {k^ - e-^)p'^{t, 1) 

+ 2{k - l)(a + kd)p{t, l)W{t) + - l)(2a + (A: - 

Let /9 > be small enough so that 

e-^ > (3.27) 

Then 

L(t) = -/3 [\-l^''p'^{t,x)dx 
Jo 



+ 



+ 



{k-l){a + kd) 



n 2 



- l)(2a + (/fc - l)d) + 



(A: - l)2(a + A;d)^ 
e-/3 - A:2 



< -/3 / e-^''p^{t,x)dx 
Jo 

+ \~ ^, , [(1 - /t)a2 + 2ad{k - e-^) + e-^(l - A;)(i21 w^2(^) 

e P — 

= -/3 C e-^''p^{t,x)dx 
Jo 



[{k - l)a + (e-^ - k)df 



Take 



e-^-k , 



which verifies ^^7I^ since jdj < 1 and ([3221) • Moreover, by ([3:23]) and (I32E1), we get 

L{t) < -I3[l-d^{e^ - [\-^'' p\t,x)dx. 



(3.28) 



(3.29) 



As /3 ^ 0+, 1 - (i2(e^ - l)2e"^^-2 ^ i _ ^2 > q_ Hence, taking /? > small enough, we may 
assume that 



l-d\e^-lfe-^r^>0 



(3.30) 



From ([3:29]) and (iOOjl . one has 



/3 



L(i) < -TT [1 - d\e^ - l)2e-^/3-2]L(t) 
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Therefore, there exists a constant a = a{p, k) > such that 

L{t) < -aL{t), yt € [0,+oo). (3.31) 
Finally, we conclude from the fact that 

m < Cs j^' e-^''pl{x)dx < C3||po|li2(o,i) (3-32) 

and p.3ip that there exists a constants C = C(p, /c) > such that p.4p holds. This concludes 
the proof of Theorem 13.11 for the case \d\ < 1. 
Case 2. d>l. 

In this case, the construction of Lyapunov function is not as direct as in the case \d\ < 1. 

Let ^ 

Vi{t):=[ p'^{t,x)dx + bW^{t), VtG[0,+oo), 
Jo 

where b S (0, +oo) is a constant to be determined. We compute the time derivative of Vi 
along the solution of ()3.3p . One has 



Vi{t) :=- [ {p\t,x)),dx + 2bW{t)W{t) 
Jo 

= p\t, 0) - p\t, 1) + 2bW{t)ip{t, 0) - pit, 1)) 

= (k^ - l)p\t, 1) + 2{k - 1)(6 + kd)p{t, l)W{t) + d{k - 1)[26 + (A; - l)d]W'^{t). 
Since k € (—1, 1), in order to then deduce that Vi{t) < 0, it suffices to require that 
4{k - l)^{b + kdf - M{k^ -l){k- l)[2b + {k - l)d] < 0, 

that is. 



Taking 
readily gives 



4{k-lf{b-df < 0. 



b := 



Let 



Vi{t)= [ p'^{t,x)dx + dW^{t), (3.33) 
Jo 

Vi{t) = {e - l){p{t, 1) + dW{t)f. (3.34) 

C{t,x) := p{t,x) + dW{t), i G (0,+(X)),x e (0, 1). (3.35) 



By (1231]) and (^M) 

Vi{t) = ik^-l)eit,l). (3.36) 
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From ()3.3p and ()3.35p . one gets that satisfies the following Cauchy problem 

' Ut,x)+^xit,x) = dW{t), t G (0,+oo),x G (0, 1), 
e(0,x) =po(x)-p + (it^(0), xG(0,l), 
yi{t,Q)=K{t,l), tG(0,+oo). 



Let 



Then 



V2{t):= I e-^'C^(t,x)dx, Vt G [0, +00) 
Jo 



(3.37) 



(3.38) 



V2{t) :=2 / e-^i{t,x)it{t,x))dx 
Jo 

= 2 [ e-^^(t, x){-^^{t, x) + dW{t))dx 
Jo 

ix=l /"l 1-1 



e-''^^{t,x)dx + 2dW{t) / e-'-'^{t,x)dx 
^=0 Jo Jo 



-e-'e{t, 1) + e{t, 0) - V2{t) + 2a!(^(t, 0) - ^(t, 1)) /' e-^CC^, x)dx 

Jo 

{k^ - e-^)f{t, 1) - V2{t) + 2d(A; - l)C{t, 1) / e-^e(t, x)dx. 

Jo 



By the Cauchy- Schwarz inequality, 

"1 2 
e-^^(t,x)dx < (1 - e-^)l/2(i)- 



From (|3.39|) and (|3.40|) . there exists a constant A = A{d, k) > such that 



V2{t)<--V2it)+Af{t,l). 



Let 



V{t) :- 



2A 



1-A:2 

Combining (pM) . (piTO and (pM) . one has 

2^ 



^i(t) + 1^2(0, VtG[0,+oo). 



V{t) 



^_j^2^l{t)+V2{t)<--V2{t), VtG [0,+oo). 



(3.39) 

(3.40) 

(3.41) 
(3.42) 

(3.43) 



Notice that 



and 



[ f{t,x)dx= [ p'^{t,x)dx + {2d + d^)W^{t) 
Jo Jo 



W'^{t) = (^j p{t,x)dx^ < j p^{t,x)dx. 

Obviously, there exist constants Bi = Bi[d) > (i = 1, 2, 3, 4) such that 

i?i||p(t,-)lli2(o,i) < B2V2{t) < V{t) < BsV2{t) < i?4||p(t,-)lli2(o,i), Vt G [0,+^)). (3.44) 

Now we conclude from (|3.43|) and (|3.44|) that there exist constants a = a{'p,k) and C = 
C(j), k) such that 

V{t) < -aV{t), yt G [0,+oo) (3.45) 
and finally ()3.4p holds. This finishes the proof of Theorem 13.11 for the case d>l. □ 
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4 Stabilization to p for the nonlinear system 

In this section, we stabilize the nonhnear system to an equihbrium p € M by using Lya- 
punov function approach. By Lemma l2. II and Lemma 12.21 it suffices to construct Lyapunov 
functions for classical solutions. We will divide our main results into two cases: p = and 
p 7^ 0. We will see later that the situation of p 7^ is much more complicated than that of 
p = which implies d = 0. 



4.1 Exponential stability of with a Lyapunov function approach 

In this subsection, we prove a stabilization result for the case that p = 0: we give explicit 
feedback laws leading to semi-global exponential stability of p = 0. 

Theorem 4.1. For every k € (—1, 1) and every R > 0, there exist constants C = C{k, R) > 
and a = a{k, R) > such that for any pQ € L^{0, 1) with 



IIpo||li(o,i) < -R) 
the solution p € C'^([0, +00); L^(0, 1)) to the Cauchy problem 

Pt{t, x) + {p{t, x)X{W{t)))^, =0, t e (0, +00), X e (0, 1), 

p{0,x) = po{x), xe{0,l), 
^u{t) = ky{t), t G (0, +00) 



(4.1) 



(4.2) 



satisfies 



pit,-)\\mo,i) <Ce "*||po||l2(o,i), VtG[0,+oo) 



(4.3) 



Proof. Still motivated by the Lyapunov functions used in [H [131 ESI [37] , we introduce the 
following Lyapunov function which is an weig hted ^2(0, 1) norm of the solution: 



L{t) := [ e-^''p^{t,x)dx, VtG[0,+oo), 
Jo 



(4.4) 



where /3 > is a positive constant to be determined. Then, along the trajectories of ()4.2 
(see also ([131) and (fT3D ). 



L 



{t)= f e-P^{p\t,x)\ 
Jo 



dx 



\{W{t)) / e-^%p\t,x)),dx 



p\{W{t))L{t) - \{W{t)) le-^'^p^t, x) 



x=l 
x=Q 



^X{W{t))L{t) + X{W{t)){p\t,0) - e-l'p'it,!)) 
mWitmt) + {X{Wm-'{k^ - e-^)y\t). 
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Since k G (—1, 1), one can choose /3 > so that e ^ > k"^ and thus 

L{t) = -pX{W{t))L{t) < 0. (4.5) 

In order to get exponential decay of the solution as t +oo, it suffices to prove the 
uniform boundedness of W{-). In fact, 

t ^ \p{t,x)\dx is a nonincreasing function. (4.6) 
Jo 

Property ()4.6p can be proved by checking that, for every r € (1, 2], 

^J'\p{t,x)rdx^O, (4.7) 
From ()4.7p . one gets that, for every r E (1,2], 



t '-^ |/3(t, x)|''(ix is a nonincreasing function, 
Jo 

which gives (|4.6p by letting r — )• 1+. From (j4.6p . one gets that 

\W{t)\ < ||po||li{o,i) < ^, Vt e [0,+oo). (4.8) 

Let us define a constant 

b := inf X(s) > 0. (4.9) 

|s|<_R 

Then we get from (gJ]), (g^D, (jM} and (gj]) that 

L{t) < -bf3L{t), Vt€[0,+oo). 

Therefore we obtain 



L{t)<me-'^'= I e-^>^(x)dxe-''^*<||po|li2(o,i) 



e 



This concludes the proof of Theorem 14.11 □ 

Remark 4.1. It is easy to see that if we let k = 0, i.e., the feedback law is chosen as 

u{t) = 0, tG(0,+oo), 

then the state reaches zero after a finite time no matter what the initial data is. This fact 
shows that zero control does drive the state to zero in finite time. 

4.2 Exponential stability of p 7^ with a Lyapunov function approach 

In this subsection, our main result is the following one. 
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Theorem 4.2. Assume that d > —1. Let k S (—1,1). Then there exist constants e = 
e{'p,k) > 0, C = C{'p,k) > and a = a{'p,k) > such that the following holds: For every 
po G i^(0, 1) with 

IIpo(-) -^IIl2(o,i) < e> (4.10) 
the weak solution p € C''([0, +oo); L^(0, 1)) to the Cauchy problem 



satisfies 



pt{t, x) + {pit, x)XiW{t)))^ =0, t£ (0, +oo), X G (0, 1), 
p{0,x) = po{x), xG(0, 1), 
[u{t) - pX{p) = k{y{t) - pX{p)), t G (0, +00) 

\\p{t, •) - p||l2(o,i) < Ce-^^Wpoi-) - p||l2(o,i), Vt G [0, +00). 



(4.11) 



(4.12) 



Proof. As we did in Section 13.11 and Section 13.21 we may assume that A(p) = 1, which, 
together with ()1.7p . yields d = 'pX'{p). 
Let 

p{t,x) := p{t,x) -p, W{t) ■.= W{t) -p, po{x) := po{x) -p, 
A(t) := A(p + W{t)), u{t) := Xit)p{t, 0), y{t) := A(t)p(t, 1). 

The system ()4.1ip is then rewritten as follows 

pt{t,x) + X{t)p^{t,x) =0, t£ (0,+oo),2; G (0,1), 

p(0,x) =po(x), xG(0,l), (4.13) 
u{t) = ky{t) + {k- l)p(X{t) - 1), t G (0, +00). 

Until the end of the proof of theorem 14. 2^ we omit the symbol~. In particular we rewrite 
(|4.13p as the following system 



pt{t, x) + X{t)p^{t, x)=0, te (0, +CX)), X G (0, 1), 
p{0,x) = po{x), xG(0, 1), 

u{t) = ky{t) + {k- l)-p{X{t) - 1), t G (0, +00), 



(4.14) 



where 



W{t)= [ p{t,x)dx, X{t) = X{p + W{t)). 
Jo 

Correspondingly, the assumption (|4.10p and the exponential decay estimate (|4.12p become 



IIpo||l2(o,i) < e> 

l|p(*r)|lL2(o,i) < (^e-^^llpollL^co,!)' ^ [0,+oo). 



(4.15) 
(4.16) 



Similar to the linear case, we divide our proof into two cases: |d| < 1 and d > 1. 
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Case 1: \d\<l. 

We define a Lyapunov function by (|3.22p . where a is given by (j3.28p and /3 > is taken 
small enough so that ()3.24p and (|3.27p hold. Then, L{t) is positive definite for every t > 
and there exist four constants Ci = Ci{d, k, f3) > (i = 1, 2, 3, 4) such that ()3.25p holds. 

Let us compute the time derivative of L{t) for any classical solution of (j4.14p . Note that 

W{t)= [ pt{t,x)dx = X{t){pit,0)-pit,l)) = u{t)-y{t). (4.17) 
Jo 

We get, from ([SJl]), ([TO]) and (glTD, that 

L{t) = -X{t) [ e-^''{p'^{t,x))^dx + 2aW{t)W{t) 
Jo 

= -/3A(t) / e-'^''p^{t,x)dx+ '^^^^^ - e V(t) ^ 2aW{t){u{t) - y{t)) (4.18) 

Jo Ht) 
1 



/3A(t) / e"'^^p^{t,x)dx + Ai, 



where 



Ai = ^'^'^ ^[^^""'^'^ + 2aW{t){u{t) - y{t)) 

_ [kyit) + ik - mXit) - l)]--e-V(t) ^ ^^^^ _ ,)^(,)[^(,) ^ _ 



A(t) 



A(t) 

(/c-l)[fcp(A(t)-l) + aA(t)]Ty(t) 

y(i) + ^^3^3^ 



{k - lfa^X\t)W^{t) 
X{t){B-e-P) 



2a{e~P -k){k- l)Xit)Wit)piX{t) - 1) + e-^{k - l)^p^(A(t) - 1)^ 

A(t)(A;2_e-/3) • ^ 

Since A is of class C^, one has 

X{t) = X{p + W{t)) = l + X'{p)Wit) + o{l)Wit), VtG[0,+oo). (4.20) 

Here and hereafter, we denote o(l) for various quantities (may be different in different situ- 
ations) satisfying the following property 

> 0, 3e > such that (^(|W"(t)| < e ^ |o(l)| < 5), Vt€[0,+oo)). (4.21) 

Then, using lK2M . (ITOD and (g^Ql), one has 

^ [{k - l)2a^ + 2a(e-^ - fc)(fc - l)d + e-^{k - lfd'^]W'^{t) + o{l)W\t) 
^- (e-/3 - A;2)(l + o(l)) 

< [d\l - e"^) + o{l)]W^{t). (4.22) 
Combining ([3:25]) . (0201), (jCT]) and (I02]l . we get 

L(i) < (-/3A(t) + d2(e/3 - l)2e-/3/3-i + o(l)) e-^'"" p^{t,x)dx 

1 
C3 



< (l - d^{e^ - ife-^r^ + 0(1) ) L{t). (4.23) 
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Prom the fact that 1 - ^^(e^ - lfe~'^l3~'^ ^ 1 - > as /? ^ 0+, we have, for /3 > 
smah enough, 

l-d^{e'^-lfe-'^p-^>0. (4.24) 

Finally, we take /3 = /3(p, A;) > small enough so that (j3.24p and (j4.24p hold. Let us also 
emphasize that, from the Cauchy-Schwarz inequality and ()3.25p . there exits C = C{'p,k) > 
such that 

\W{t)\^<CL{t), (4.25) 

which, together with (|4.15p . ()4.2ip . ()4.23p and ()4.24p concludes the proof of Theorem Ol in 
the case \d\ < 1. 
Case 2: d>l. 

In this case, we are going to prove that the Lyapunov function in the type of ()3.42p still 
works for the nonlinear control system ()4.14p . 

Let's first compute the time derivative of Vi{t) (see (|3.33p for definition), using ()3.35p . 
(mZl) and 

Vi{t) = -X{t) [ {p^{t,x))^dx + 2dW{t)W{t) 



= X{t)[p\t, 0) - p\t, 1) + 2dW{t){p{t,0) - p{t, 1))] 

- i^y(^)^i^-mm-i)?-yHt) ^ ^^^^ ,^wmyit)^-pm - d) 

= X{t){k^ - l)(p(t, 1) + dW{t)f + [o{l)W{tMt, 1) + o{l)W\t) 

= x{t){k^ - i)e{t, 1) + oii)w{t)at, 1) + o{i)w\t). 

Thanks to the Cauchy-Schwarz inequality and (|3.44p . it follows that 

Vi<{k^-l + 0(1)) fit, 1) + o{l)W\t). (4.26) 
By the definition of ^, it is easy to get that ^ satisfies the following Cauchy problem 
Ct{t, x) + X{t)C^{t, x) = dW{t), t G (0, +CX)), x G (0, 1), 

C{0,x) = po{x) + dW{0), xG(0,l), (4.27) 
e(t, 0) = kat, l) + ik- 1) (^^^f^ - dW{t)^ , t G (0, +oo). 

Let, again, V2 be defined by (|3.38p . Then the time derivative of V2 satisfies: 



V^{t) = -X{t) [\-''{f{t,x))^dx + 2dW{t) C e-^'iit, 
Jo Jo 

= -X{t) [ e-''f{t,x)dx - X{t)\e^^f{t,x)Y^^ 



x)dx 







+ 2(iA(t)(e(t,0) -e(t,l)) / e-''^{t,x)dx. 

Jo 
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Hence 

Mt) = -x{t)V2it) - x{t) 



e-'e{t, 1) - ( kC{t, l) + {k- 1) (^^^^^ - dW{t) 



+ 2d{k - l)\{t) 1) + ^^^^^^^^^ - dW{t)^ e-^'e(t, x)dx 

= -xm2it) + A(t)(fc2 - e-')eit, 1) + oii)w{t)at, i) + oii)w\t) 

+ 2dik-l)\{t){C{t,l) + o{l)W{t)) [ e-''C{t,x)dx. (4.28) 

Jo 

By the Cauchy-Schwarz inequality, (|3.40p . (|3.44p and (|4.28p . one has, for A = A{p,k) > 
sufficiently large, 

Mt) < (-1 + o{l))V2{t) + ^(1 + o{l))e{t, 1) + o{l)W\t). (4.29) 
We stiU define V by ([332]). Combining ()i:26]) and (|09]l . one has 
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(4.30) 



< (-1 + o{l))V2{t) + {-A + o(l))e'(t, 1) + o(l)i^2^t). 
Consequently, there exists a = a{'p,k) > such that 

Vit) <{-a + oil))V{t). (4.31) 

This, together with ([02]) ([OI]) . (fiJ5]) and dMH), concludes the proof of Theorem K2\ in 
the case d > 1. □ 
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